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Nuclear matrix elements (NME) for exchange of arbitrary mass neutrinos are calculated in the 
interacting boson model (IBM-2). By combining the NME with the phase space factors (PSF), 
expected half-lives for neutrinos of mass mM and coupling UeN are estimated. Limits on sterile 
neutrinos with masses in the eV, keV, MeV-GeV, and TeV range are given. 
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I. INTRODUCTION 

In recent years, the possible occurrence of sterile neu¬ 
trinos with mass rriN in addition to the three known 
species mi, m 2 , m 3 , has attracted considerable attention, 
and searches are under way to detect their presence in 
oscillation experiments and in accelerator experiments. 
Sterile neutrinos were introduced by Pontecorvo jl| in 
1968 as neutrinos with no standard model interaction. 
Several suggestions have been made for, among others, 
sterile neutrinos in the eV mass range [2, Qj.,™ the keV 
mass range , in the MeV-GeV mass range p, , and in 
the TeV mass range Q- 

Sterile neutrinos, if they exist, will contribute to neu¬ 
trinoless double beta decay (DBD). It is therefore of in¬ 
terest at the present time to estimate the expected half- 
life for Majorana neutrinos of arbitrary mass. In a previ¬ 
ous series of papers we have considered the case of very 
light, rriN pf, and very heavy misi ^ pf, neutrinos, 
where pF lOOMeV is the nucleon Fermi momentum in 
the nucleus (in this article we use units c = 1). In these 
cases the half-life factorizes to 

[^ 1 / 2 ]-' = Go. \Mo.\^ |/(m„C/eOl' , ( 1 ) 

where Gqi, is a phase space factor (PSF), Mq^, the nu¬ 
clear matrix element (NME), and / is equal to, for light 
neutrino exchange, 

/ = {m,,) = ^ {Uekf mk, ( 2 ) 

rup 

k—light 

and, for heavy neutrino exchange 

/ = mp(m;J), (m-i)= X 

I L '^kh 

kh—heavy 

( 3 ) 
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When the mass mjv is intermediate, and especially, when 
it is of the order of magnitude ofpF, the factorization ( 1 ) 
is not possible, and physics beyond the standard model is 
entangled with nuclear physics. In this case, the half-life 
can be written as 

= Go. 


The matrix elements Mo.(mAr)mAr/me have been cal¬ 
culated within the framework of the interacting shell 
model, ISM, @ and quasiparticle random phase approx¬ 
imation, QRPA, Here we present results within the 
microscopic interacting boson model, IBM-2, with isospin 
restoration From the NME and the PSF that we 
have previously tabulated ini , we then make estimates 
for half-lives and set limits on some of the suggested ster¬ 
ile neutrino species. 


^' TOr 


N 


( 4 ) 


II. FORMALISM 


Although several formulations of Or^/3/3 decay are avail¬ 
able, we use here the formulation of Simkovic et al. fl^ , 
as adapted to the case of neutrinos of arbitrary mass, 
ruN i- The transition operator T{p) is the same as in 
Eq. (5) of our previous paper with f = rriFi/me and 


v{p) 


2 1 
^ + m|f -b AJ 


( 5 ) 


where v{p) is the so called neutrino "potential". The 
product fv(j>), 


fy{p) 


TTlAT 2 1 

(^y/p"^ + m% -b 


( 6 ) 
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has the limits: 


with 


rriM —t 0 


rriN —t oo 


fv 


fv 


rriN 2 1 

TOe TT p 

rrii^ 2 1 TTip / 2 1 

TOe TT 'TTIn \TT nielTlp 


(7) 


as in our previous calculations for light and heavy neutri¬ 
nos [l^ . In Eq. (5) A is the closure energy assumed 
to be a smoothly varying function of the mass number 
A as given in Table 25 of Ref. H. All other quan¬ 
tities, form factors h(p), momentum dependence of the 
cou plin g constants gv(p^) and are the same as in 

floL Also, in this paper we use Argonne short-range 
correlation fl5| as in our previous article 0- 

The nuclear matrix elements appearing in Eq. (4) can 
be written as 

Mo„{mN) imN), 

{rriN) =mXX {’’tT'n) - MX''\mN) ( 8 ) 


with the axial vector coupling constant gA explicitly writ¬ 
ten in front, where is the dimensionless ma¬ 

trix element calculated with the neutrino potential of 
Eq.(5). 


III. RESULTS 


In Table |T] we show the dimensionless IBM-2 nuclear 
matrix elements { 171 ^) 171 ^/nie for Oz//3/3 decay with 
Argonne-SRC and gv/gA = 1/1.269 to the ground state, 
0 ^, and first excited 0^ state, Oj, broken down into 
Gamow-Teller, GT, Fermi, F, and Tensor, T, contribu¬ 
tion and their sum, according to Eq. (8) for several val¬ 
ues of rriN- These matrix elements are plotted in Fig. 
[T] for three isotopes, ^®Ge, ^^°Te, and ^^®Xe of interest 
to GERDA, GUORE and KamLAND-Zen/EXO experi¬ 
ments. 

The interesting aspect of Fig. [T] is that the curves 
peak at mAT ^ lOOMeV, the scale set by the nucleon 
Fermi momentum in the nucleus, pp. If sterile neutrinos 
of this mass exist, their contribution to the half-life is 
enhanced. The relatively simple behavior of Fig. [T]allows 
one to write a simple formula to estimate the effect of 
sterile neutrinos on half-lives for 0i//?/3 decay. The simple 
formula is 



Go.g\ 




mp^iUeNf , 2 ’ 

(9) 


where and are the matrix elements calcu¬ 

lated in the limits mpf —> 0 and ttln ^ 00 , respectively. 
The simple formula ® is identical to that given in @ 
and suggested in fl^ . It provides a good approximation 
to the exact calculation, except in the region mpi pf, 
as one can see in Fig. [T| where the simple formula dH]) 
is compared with the exact calculation. The values of 
and were given in our previous calculation 

M and are reproduced in Table El together with the 
value of for all nuclei of interest. The error in this 
table has been increased from 16% and 28% in 0 to 
32% and 56% for and respectively, to ac¬ 

count for a SRG different from Argonne. SR C g reatly 
affect the value of as discussed in Ref. |l3l |. 

It is of interest to compare our results with those of 
the interacting shell model (ISM) jH and the quasiparti¬ 
cle random phase approximation (QRPA) Q. In Fig. m 
we show this comparison. All three calculations have the 
same behavior. However, the maximum of the curves, 
related to the value of the Fermi momentum, pf, is dif¬ 
ferent. For example, in ^®Ge the maximum of the IBM-2 
curve is at ^ lOOMeV, while for QRPA-Tii it is at 
160MeV and for ISM at 180MeV. While, for tun < Pf 
QRPA and IBM-2 agree, for mN > Pf the agreement is 
between IBM-2 and ISM. The actual form of the curve 
arises from the interplay between short range correla¬ 
tions, SRC, and the wave functions close to the origin. 
Also, while the tensor component Mp plays a small role 
for TnN < Pf, it plays a large role for mAr > pp, as one 
can see from Table I. Because of the divergence of the 
tensor component at the origin, the numerical evaluation 
of its matrix elements is very sensitive. For this reason 
we have increased the IBM-2 error estimate of in 

Table II. 


IV. HALF-LIVES AND LIMITS ON STERILE 
NEUTRINO CONTRIBUTIONS 


Using Eq. (4) and our previously calculated we 
can calculate the expected half-life for a single neutrino 
of mass mAT with coupling UeN- This is shown in Fig. [3] 
for several values = 10“^ —10“®. Note the resonant¬ 
like behavior at to at ~ Pf- For other values the 
half-lives scale as In the figure, the current lim¬ 

its from GERDA 0 in ™Ge, CUORE-0 0 in ^^^Te, 
KamLAND-Zen |l9l |. and EXO in ^^®Xe are also 
shown. The situation is summarized in the exclusion plot 
di of Fig. m 

The calculated half-lives in Fig. [3] are for gA = 1.269. 
For other values of gA they scale as f/gX The renormal¬ 
ization of the axial vector coupling constant, gA, in nuclei 
is, at the present time, a major issue. Three possibilities 
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TABLE 1. IBM-2 nuclear matrix elements (dimensionless) for Qvfi P decay with Argonne SRC, QvIqa = 

1/1.269, as a function of mjv in ^®Ge, ^^°Te, and ^®®Xe. 


A 

rriN 

(GeV) 



Of 




Of 


Lrl 

rriN 

7iy(0!^) rriN 

miL 

Lrl 

7i/r(0'^) rriN 

rnN 


^^Ge 

0.00001 

0.0891 

-0.0134 

-0.00494 

0.0925 

0.0386 

-0.00520 

-0.00203 

0.03980 

76Ge 

0.0001 

0.891 

-0.134 

-0.0494 

0.92 

0.386 

-0.0520 

-0.0203 

0.398 

"^Ge 

0.001 

8.89 

-1.34 

-0.494 

9.22 

3.58 

-0.520 

-0.203 

3.97 

^^Ge 

0.010 

81.4 

-13.3 

-4.92 

82.8 

35.1 

-5.16 

-2.01 

36.3 

^^Ge 

0.100 

320 

-91.0 

-36.4 

340 

130 

-34.5 

-14.3 

137 

^^Ge 

1.000 

108 

-39.5 

-26.6 

106 

40.1 

-13.8 

-9.66 

39.1 

^^Ge 

10.00 

9.59 

-4.01 

-2.95 

9.14 

3.53 

-1.40 

-1.07 

3.34 

^^Ge 

100.0 

0.957 

-0.402 

-0.295 

0.911 

0.353 

-0.140 

-0.107 

0.333 

76Ge 

1000 

0.0957 

-0.0402 

-0.0295 

0.0911 

0.0353 

-0.0140 

-0.0107 

0.0333 

i30Te 

0.00001 

0.0681 

-0.0127 

-0.00285 

0.0731 

0.0499 

-0.00878 

-0.00174 

0.0536 

i30Te 

0.0001 

0.681 

-0.127 

-0.0285 

0.73 1 

0.499 

-0.0878 

-0.0174 

0.536 

i30Te 

0.001 

6.79 

-1.27 

-0.285 

7.30 

4.97 

-0.878 

-0.174 

5.34 

i30xe 

0.010 

62.6 

-12.6 

-2.83 

67.6 

45.6 

-8.70 

-1.73 

49.2 

i30Te 

0.100 

253 

-82.8 

-21.2 

284 

172 

-54.9 

-12.5 

194 

i30Te 

1.000 

88.0 

-34.9 

-16.4 

93.3 

53.3 

-21.0 

-9.18 

57.2 

i30Te 

10.00 

7.81 

-3.55 

-1.82 

8.19 

4.70 

-2.12 

-1.02 

5.00 

i30xe 

100.0 

0.779 

-0.356 

-0.183 

0.818 

0.469 

-0.212 

-0.102 

0.499 

i30xe 

1000 

0.0779 

-0.0356 

-0.0183 

0.0818 

0.0469 

-0.0212 

-0.0102 

0.0499 

i36Xe 

0.00001 

0.0562 

-0.0102 

-0.00216 

0.0603 

0.0294 

-0.00468 

-0.000755 

0.0316 

i36Xe 

0.0001 

0.562 

-0.102 

-0.0216 

0.603 

0.294 

-0.0468 

-0.00755 

0.316 

i36Xe 

0.001 

5.60 

-1.02 

-0.216 

6.02 

2.94 

-0.468 

-0.0755 

3.15 

i36Xe 

0.010 

51.5 

-10.1 

-2.16 

55.6 

26.7 

-4.64 

-0.755 

28.8 

i36Xe 

0.100 

204 

-65.8 

-16.1 

229 

93.4 

-28.6 

-5.43 

106 

i36Xe 

1.000 

69.4 

-27.4 

-12.5 

73.8 

26.4 

-10.3 

-4.02 

28.8 

i36Xe 

10.00 

6.15 

-2.79 

-1.40 

6.48 

2.32 

-1.03 

-0.448 

2.51 

i36Xe 

100.0 

0.614 

-0.279 

-0.140 

0.647 

0.232 

-0.103 

-0.0449 

0.251 

i36Xe 

1000 

0.0614 

-0.0279 

-0.0140 

0.0647 

0.0232 

-0.0103 

-0.00449 

0.0251 





FIG. 1. (Color online) IBM-2 NMEs for neutrinos of arbitrary mass plotted as a function of niN in a) ^®Ge, b) ^^°Te, and c) 
^®®Xe. Blue squares represent the exact calculation for niN = O.OOlGeV, O.OlGeV, O.lGeV, IGeV, lOGeV, joined together by a 
Mathematica interpolating formula. The curve is obtained using the simple formula ®. 


are often considered, gA = 1.269 (free valne), gA = ^ ing). Accordingly, the excluded region varies in Fig. O 
(qnark valne), and gA = 1.269A“*^'^® (maximal qnench- From this fignre, one can see that if gA is renormalized 
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FIG. 2. (Color online) Comparison between IBM-2 (Argonne SRC) (red), QRPA-Tii (Argonne SRC) 3 (blue), and ISM 
(UCOM SRC) (green) NMEs for neutrinos of arbitrary mass plotted as a function of mjv in a) ^®Ge, b) ^®°Te, and c) ^®®Xe. 



FIG. 3. (Color online) Expected half-life for a single neutrino of mass rriN with coupling = 10“^ — 10“® and qa = 1.269 
for a) ^®Ge, b) ^®°Te, and c) ^®®Xe. Blue squares represent the exact calculation for mjv = O.OOlGeV, O.OlGeV, O.lGeV, IGeV, 
lOGeV. The smooth curve is obtained using the simple formula ([9|. The experimental limits from CERDA [3, CUORE-0 
[l^ . KamLAND-Zen fl^ . and EXO are also shown. The excluded zone is that below these limits. 



mA?(GeV) 


FIG. 4. (Color online)Excluded values of |t/ejv|^ and mjv in 
the mjv-|t/ejv|^ plane, for qa = 1.269. For each experiment, 
CERDA [I3|, CUORE-0 [l||, KamLAND-Zen 0, and EXO 
[20II . a band of values is given, corresponding to our error 
estimate. 



mA?(GeV) 


FIG. 5. (Color online)Excluded region for different values of 
Qa, for ^®®Xe (EXO collaboration f^ l. 


in OvjSP as much as in 2i//3/3 (Til, fl^. and as in single-/? 
decay the limits in Fig. [3] should be multiplied 

by factors of ^ 22 in ^®Ge, of 33 in ^^°Te, and of 34 
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TABLE 11. The matrix elements and the value 

of (p^) for all nuclei of interest in DBD. 


Decay 


jVf(oi^h) 

(MeV) 

48Ga 

1.75(56) 

47(26) 

113 

"^Ge 

4.68(150) 

104(58) 

103 

82Se 

3.73(120) 

83(46) 

103 

962r 

2.83(90) 

99^5) 

130 

looMo 

4.22(135) 

164(92) 

137 

llOpd 

4.05(130) 

154(86) 

135 

“6Cd 

3.10^00) 

11072 ) 

130 

i24Sn 

3.19(102) 

79(44) 

109 

128Xg 

4 . 10731 ) 

101(57) 

109 

130Xg 

3.707i8) 

92(52) 

109 

I34xe 

4.05^30) 

91(51) 

104 

i36Xe 

3.05(98) 

73(41) 

107 

148Nd 

2.31(74) 

103(58) 

146 

150Nd 

2.67^5) 

11675) 

144 

i^^Sm 

2.82(90) 

11373 ) 

139 

leoGd 

4.08(131) 

15577 ) 

135 

198pt 

2.19(70) 

10478 ) 

151 

232Th 

4.04(128) 

15979 ) 

138 

238^ 

4.81^54) 

189(106) 

137 


in i36Xe. 


V. LIMITS ON SUGGESTED STERILE 
NEUTRINOS 


As mentioned in the introduction several types of ster¬ 
ile neutrinos have been suggested. We consider here a 
family of neutrinos at the eV scale ii, a family of neu¬ 
trinos at the keV scale Q, one at the MeV-GeV scale 
d, @ , and one at the TeV scale 0 . 

The total contribution to the half-life can then be 
rewritten, using Eq. m 


[^ 72 ]-^ = Go.gl 


■ Glk^rik + —Y1 ^eirrii + 

• ^ m _ < ^ m _ t ^ 


fc=l 




■T.'Jl 


rriN 


P Z_^ '^eN _|_ ^2 


N 


N 


^ 1 
y] uL — 


kh = l 




( 11 ) 


Here we have separated the contribution of the light, 
rriN Pf, neutrinos, into known k = 1,2,3, unknown 
at eV scale, i, unknown at keV scale, j, and used the 
expression appropriate for them in terms of _ We 

have also explicitly written the contribution of intermedi¬ 
ate mass, mjv ~ Pf, neutrinos at MeV-GeV scale in the 
form Q, and finally added the contribution of heavy, 
mj\! ^ Pf, neutrinos at the TeV scale, using the form 
appropriate for them in terms of 


To set limits on sterile neutrino contribution, one 
should do a simultaneous analysis where all contribu¬ 
tions are included. However, in view of the structure 
of Eq. (11), it is sufficient to consider, for the purposes of 
this paper, separate contributions of light, intermediate, 
and heavy neutrinos. This separation may, in general, 
not be possible, as the separate contributions may inter¬ 
fere with each other. 


The analysis of light (m^v Pf) sterile neutrinos can 
be done as in previous publications [l^ . Introducing 


the quantity 

3 

{mN,Hght) = Glk^k + y] Ul^m, + Y2 Uljrrij (12) 

k—1 i j 


and comparing the experimental limits with the calcu¬ 
lated half-lives using 



G 0.4 


/ {jTIN, light) 

V kUe 




(13) 


we obtain the results in Table Hill Gonsidering, for exam¬ 
ple, the case suggested in 0 of a 4th neutrino with mass 
1714 = leV and \Ue 4 \ = 0.03, we have 


3 

{rriN,light) = ^ (14) 

k^l 

where the unknown phase is 0 < 04 < 27r. Using for the 
contribution of the known neutrinos the best fit values 
given in j^, we obtain the result in Fig. [ 6 l The effect of 
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the 4th neutrino is to add to the average mass a contribu¬ 
tion of 30meV, making the spread of the allowed values 
in Fig. [6] larger than in the corresponding Fig. 7 of 0 
without 4th neutrino and thus improving the possibility 
to detect it in the next generation experiments. 



FIG. 6. (Color online) Current limits for {mM.uaht) from 
the CUORE-0 [III, IGEX [U, NEMO-3 [ill, KamLAND- 
Zen [Hi, EXO lloll . and CERDA [ItII experiments, and IBM-2 
Argonne SRC NMEs and qa = 1.269. The value of Ref. [H 
is shown by X. The figure is in semilogarithmic scale. Red 
shows the normal hierarchy and green the inverted hierarchy. 
In this figure the scenario suggested in [3|, relevant to LSND 
and reactor anomaly, is considered. 


The combined analysis of heavy, rriN 'S> pf, and inter- 
mediate, my pf, mass neutrinos can also be done as 
before (ift [l3l | . The inverse half-life is given by 



Go.g\ 




br 


(15) 


with 


TABLE III. Left: Calculated half-lives, Eq. (13), in IBM- 
2 Argonne SRC, gA = 1.269, and = leV. Right: 

Upper limit on sterile neutrino mass {mM,light) from current 
experimental limit from a compilation of BarabcishjM. The 
value reported by Klapdor-Kleingrothaus et al. [H, IGEX 
collaboration [H, and the recent limits from KamLAND-Zen 
[Hi EXO [Hi and CERDA fl^ are also included. 


Decay 

rO;,(1024yr) 

U°/ 2 .exp(yi') {m-N, light) (eV 

4®Ca^-4®Ti 

1.33 

> 5.8 X 10^2 

< 4.8 

76Ge^76se 

1.95 

> 1.9 X 1028 

< 0.32 



1.2 X 1028_^ 

0.40 



> 1.6 X 1028_^ 

< 0.35 



> 2.1 X 1028'= 

< 0.30 

82Se-^®2Kr 

0.71 

> 3.6 X 1023 

< 1.4 

96Zr^96Mo 

0.61 

> 9.2 X 1021 

< 8.1 


0.36 

> 1.1 X 1024 

< 0.57 

iiopd_^iiOGd 

1.27 




0.63 

> 1.7 X 1023 

< 1.9 

124Sn^l24Tg 

1.09 



i28Te-)>428Xe 

10.19 

> 1.5 X 1024 

< 2.6 

i30Te-^i30xe 

0.52 

> 2.8 X 1024 

< 0.43 


10.23 




0.74 

> 1.9 X 1028_f_ 

< 0.20 



> 1.1 X 1028® 

< 0.25 

148Nd-).148Sin 

1.87 



i50Nd-).i50sni 

0.22 

> 1.8 X 1022 

< 3.5 


4.19 



i™Gd^i®°Dy 

0.63 



198pt_^198jjg 

2.77 



232Th-^232u 

0.44 



238U^238p^ 

0.13 




=>• Ref. 
Ref. 
Ref. 
a Ref. 
® Ref. 



maximal quenching, they should be multiplied by larger 
a factor which depends on A. For ^®Ge, this factor is 


~4.7. 


77 = = nip 


- 


rrik^ 


E 

N 


u. 


rriN 


eN 


(p2) 




( 16 ) 

By comparing the experimental limits with the calcu¬ 
lated half-lives, we obtain the results in Table IV. The 
effect of additional intermediate mass neutrinos is to add 
a contribution to tj. 

All the above limits have been set with gA = 1.269. 
As mentioned in the previous section, gA is renormalized 
in models of nuclei • The value of gA appears 

to the fourth power in the half-life. If = 1, the 

quark value, all limits for {mFi^ught) in Table III should be 
multiplied by a factor of ^ 1.6. 115 ^^ 6 ^/”^ = 1.269A“°-^®, 


VI. CONCLUSIONS 

In this article, we have calculated nuclear matrix 
elements for exchange of neutrinos of arbitrary mass 
within the framework of the microscopic interacting bo¬ 
son model-2 in all nuclei of interest to double beta decay. 
A simple analytical formula has been shown to de¬ 
scribe the full calculation well. Using this formula and 
the experimental limits on half-lives, we have then set 
limits on possible sterile neutrino contributions in neu¬ 
trinoless double beta decay. As a result, a good fraction 
of the parameter space in the mN-\UeN\ plane is ex¬ 
cluded, especially for neutrinos in the intermediate mass 
range, Figm 
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TABLE IV. Same as Table Hill for Eq. (16) with 77 = 1 x 10 The last column shows the lower limit on {mM,heavy)- 


Decay 

U°/2 (lO^Vr) 


1,71(10-3) 

{tUN, heavy) (GcV) 


0.72 

> 5.8 X 10^2 

<0.36 

> 2.6 

^Ge^-^^Se 

1.51 

> 1.9 X 1025 

<0.028 

>33.5 



1.2 X 1025^ 

0.035 

26.8 



> 1.6 X 1025^ 

<0.031 

>30.3 



> 2.1 X 1025_^ 

< 0.027 

>34.8 

®2Se-i>®2Kr 

0.55 

> 3.6 X 1023 

<0.12 

>7.82 

96Zr^96Mo 

0.19 

> 9.2 X 1021 

<0.46 

>2.04 


0.09 

> 1.1 X 1024 

<0.028 

>33.5 

iiopd_^iiOGd 

0.33 





0.19 

> 1.7 X 1023 

<0.11 

>8.53 

i24sn^i24Te 

0.67 




128pg^l28Xe 

6.43 

> 1.5 X 1024 

<0.21 

>4.47 

i30pg^i30xe 

0.32 

> 2.8 X 1024 

<0.034 

>27.6 


7.73 





0.50 

> 1.9 X 1025^ 

<0.016 

>58.6 



> 1.1 X 1025« 

<0.021 

>44.7 

148Nd^l48Sm 

0.36 




i50Nd^i50sm 

0.05 

> 1.8 X 1022 

<0.16 

>5.9 

154Sm^l54Qjj 

1.00 




i60Gd^i60p)y 

0.17 




198pt^l98Hg 

0.48 




232Th_>,232-(j 

0.11 




238U^238p^j 

0.03 





Ref. 
Ref. 
<= Ref. 
d Ref. 
Ref. 
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